o 

>v 



x 



First result from the Double Chooz reactor-neutrino experiment 

T. Matsubara on behalf of the Double Chooz collaboration 
Department of Physics, Tokyo Metropolitan University, Japan 



f"^ , We report first results of a search for the non-zero neutrino mixing angle #13 from the Double 

■ Chooz experiment. Double Chooz aims to measure the mixing angle based on i> e disappear- 

ance as a consequence of neutrino oscillation. A new generation of u e detector having f m 3 
fiducial volume is located 1 km from the two 4.25 GWth reactors at the Chooz Power Plant 
in France. Physics data taking has been continuing since April 201 1. A ratio of observed-to- 
predicted event rate of 0.944 ± 0.016 (stat) ± 0.040 (syst) was obtained in 101 days of detector 
C— running. Analyzing both the rate and their energy spectral shape, we found sin 2 2#i3 = 0.086 

± 0.041 (stat) ± 0.030 (syst) at Am 2 atm = 2.4 x 10" 3 eV 2 . 

1 Introduction 

> 

. Prior to one year ago, the neutrino mixing angle #13 had been an unknown parameter in the 

neutrino oscillation framework. The best constraint on #13 was come from the CHOOZ reactor- 
neutrino experiment^, sin 2 26>i 3 < 0.15 (90% C.L. at Am 2 atm = 2.5 x 10~ 3 eV 2 ). The T2K and 
MINOS accelerator experiments recently reported a sign of non-zero #13 with v e appearance^^. 
First reactor experiment to present a result was Double Chooz, which reported an indication 
for a disappearance of reactor v e in November 2011^. Other two reactor-neutrino experiments, 
Daya Bay and RENO, gave consistent results with higher sensitivity in early 2012^^. These 
results supporting relatively large #13 value, sin 2 2#i3 ~ 0.1, potentially encourage us to measure 
a leptonic CP violation in near future experiments. This paper describes the first results of a 
search for the non-zero neutrino mixing angle #13 from the Double Chooz experiment. 

Neutrino oscillation occurs as a consequence of non-zero mass and a mixing of mass eigen- 
states and flavor eigenstates. In the reactor-neutrino experiments, survival probability of v e is 
given by: 

P surv ~ 1 - sin 2 26 13 sin 2 (1.27Am 2 tm L/£), (1) 

where Am^ tm (eV 2 ) is the atmospheric squared mass difference, which is precisely measured by 
the MINOS experiment^. L is the reactor-to-detector distance in meters and E the v e en- 
ergy in MeV. This formula indicates the sin 2 2#i3 measurability by observing v e disappearance 
at appropriate baseline. We therefore placed a v e detector located at ~1050m under a ~300 
m.w.e. rock overburden from the two 4.25 GWth reactors at the Chooz Power Plant in France. 
In order to reduce systematic uncertainties such as neutrino flux, we will operate an identical 
detector located at ~400 m baseline under a ~120 m.w.e. rock overburden by 2013. The P e 's are 
detected through the inverse beta decay (IBD) reaction: v e + p — > e + + n. Detector based on 
hydrocarbon liquid scintillator provides free proton target and reacts with v e . Positron ioniza- 
tion and annihilation (l~8MeV) then creates a prompt signal. Neutron capture on Gadolinium 



(~8 MeV) creates a delayed signal. The signature of IBD reaction is identified by a time coin- 
cidence of r ~ 30 [is between those signals. In this reaction, the v e energy can be reconstructed 
with prompt energy as: E prompt = E(kin.) e + + 2m e ~ E Pe - (M n - M p ) + m e ~ E Pe - 0.782 MeV. 

A new generation of v e detector for the Double Chooz experiment (Figure [TJ) consists of a 
main detector, an outer veto and calibration devices. The main detector is separated into four 
concentric cylindrical tanks. Innermost 8 mm thick transparent acrylic vessel is called neutrino- 
target region. The region is filled with a liquid scintillator with a mixture of n-dodecane, PXE, 
PPO, bis-MSB and 1 g/1 Gadolinium as a beta-diketonate complex. This structure contains a 
fiducial volume for the neutrino events within 10 m 3 of the target vessel. The composition is 
chosen for radiopurity and long-term stability. Gamma-catcher region surrounds the target vessel 
to detect 7-rays escaped from neutrino-target region. It consists of 12 mm thick transparent 
acrylic vessel containing 22.3 m 3 of liquid scintillator. Buffer region is located outside of the 
gamma-catcher region. The region is filled with 110 m 3 of mineral oil to shield 7-rays from the 
Photo Multiplier Tubes (PMT) and rocks. Total 390 10-inch low-radioactive PMTs are equipped 
at the inner wall of stainless steel buffer tank to observe light from the inner volumes. These 
three layers constitute the inner detector (ID). The ID is surrounded by the optically separated 
region of the inner veto (IV) equipped with 78 8-inch PMTs. 90 m 3 liquid scintillator is filled in 
a steel tank for muon veto and shielding of spallation neutrons produced outside of the detector. 
The detector is surrounded by 150 mm of demagnetized steel to reduce a contamination due to 
external 7-rays. Outer veto system lies on main detector, which is not used for this analysis. 
The main detector is read out with 500 MHz flash- ADC electronics with customized firmware 
and a deadtime-free acquisition system. 
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Figure 1: A schematic view of the Double Chooz detector. 



2 Neutrino prediction and data analysis 



The analysis is based on 101.5 days of data including 16 days with one reactor off and one day 
with both reactors off. The expected number of u e 's are predicted as: 

where N p is the number of protons in the detector, e the detection efficiency, Pth(t) the thermal 
power of reactor, (Ef) the mean energy per fission and (af) the mean cross-section per fission. 
Mean energy per fission and associated errors are evaluated with two reactor simulation codes, 
MURE 8 and DRAGON®. Mean cross-section per fission and associated errors are estimated 
with updated reference spectra^! Neutrino prediction is based on Bugey4 measurement^^ as 
an anchor point of the mean cross-section with correction to Chooz reactor. Total systematic 
uncertainty on reactor is estimated to be 1.8%. 

Energy measurements for data analysis are based on the total charge collected by the 
PMTs. The energy is reconstructed scaling the total charge by a constant corresponding to 
~200 p.e./MeV, which is adjusted with the 2.2 MeV energy peak of neutron capture on H at 
the target center. Our Monte Carlo simulation based on GEANT4 is used to calculate the de- 
tector response. Comparison between actual and simulated calibration data has been done to 
correct the simulation and estimate an associated uncertainty to be 1.7%, using two parametric 
functions with respect to energy and position. 

We applied the following criteria to select v e candidates. Triggers within a 1000 /is after a 
cosmic muon crossing the ID or IV are vetoed to suppress spallation neutrons and cosmogenic 
backgrounds. This requirement is followed by five selections: (1) a cut rejecting events caused 
by sporadically glowing PMT bases, resulting in light localized to a few PMTs and spread out in 
time: Q ma x/Qtot < 0.09 (0.06) for the prompt (delayed) energy and rms(t S i ar t) < 40 ns, where 
Qmax is the maximum charge recorded by a single PMT and rms(t S i ar i) is the standard deviation 
of the times of the first pulse on each PMT; (2) 0.7 MeV < E prompt < 12.2 MeV; (3) 6.0 MeV 

< ^delayed < 12.0 MeV; (4) 2 [is < Ai e + n < 100 /xs; (5) a multiplicity cut to reject correlated 
backgrounds defined as no additional valid trigger from 100 fis preceding the prompt candidate 
to 400 fj,s after it. In total, 4121 neutrino candidates were observed, which is equal to 42.6 ± 0.7 
events/day on average (Figure [2]). 

Backgrounds mimicking the time coincidence have been estimated as follows. Accidental 
background caused by uncorrelated coincidence, for example, prompt event from radioactivity 
and delayed event from cosmogenic neutron capture. This background rate is estimated by 
sequentially shifted off-time window, leading 0.33 ± 0.03 events/day. Fast neutron induced by 
muon traveling the rock can interact in the detector producing recoil proton. Thermalized and 
captured neutron accompanying the recoil proton then mimics the IBD events. We estimated 
the rate to be 0.83 ± 0.38 event/day by modifying the E prompi selection criteria to be 12.2 MeV 

< ^prompt < 30 MeV and extrapolating to the signal region assuming flat energy spectrum. 
Spallation product of 9 Li induced by energetic muons emits n and /3, simulating IBD event. 
This background is studied by time distribution between a muon energy deposition > 600 MeV 
and neutrino candidates. Fitting the time distribution by a flat component and an exponential 
with the 9 Li lifetime results in 2.3 ±1.2 events/day. We account for a shape uncertainty between 
some decay branches. 

We had a chance to take data with both reactors off for ~24 hours. It was an unique op- 
portunity to validate our background estimation. Two candidates in neutrino energy window 
following high energy muon were observed, which are compatible with candidates of 9 Li back- 
ground. This result is consistent with the estimated number of BG events (Total 3.46 ± 1.26 
events / day) . 




Figure 2: Observed and expected neutrino event rate. 



3 Oscillation analysis 



Oscillation analyses have been done based on x 2 estimator. Uncertainties used for the oscil- 
lation analyses are summarized in Table [TJ A ratio of observed-to-predicted events of 0.944 
± 0.016 (stat) ± 0.040 (syst) was observed, corresponding to sin 2 26»i 3 = 0.104 ± 0.030 (stat) 
± 0.076 (syst) at Am 2 3 = 2.4 x 10~ 3 eV 2 . The analysis can be improved by using spectral in- 
formation. The positron spectrum is divided into 18 variably sized energy bins from 0.7 MeV 
to 12.2 MeV. In order to introduce bin-to-bin correlations, we use four covariance matrices to 
include uncertainties on statistics, reactor, detector and background spectral shape. Shape dif- 
ference between signal and background appears in the region from 8 to 12.2 MeV, reducing the 
uncertainties due to correlated backgrounds. Analyzing both the rate and their energy spectral 
shape, we found sin 2 26»i 3 = 0.086 ± 0.041 (stat) ± 0.030 (syst) with xV n -d-f- = 23.7/17. Ob- 
served and predicted positron energy spectra for no oscillation and the best-fit sin 2 2#i3 including 
background are shown in Figure El No oscillation hypothesis is ruled out at the 94.6 %C.L., 
which can be interpreted as an indication of the non-zero #13. This result is compatible with 
follow-up results from the Daya Bay and RENO experiments. A combined analysis with the 
T2K and MINOS accelerator experiments on #13 and CP violation phase 5 vs #13 plane for 
normal mass ordering is presented as shown in Figure 01 



4 Conclusion 



Double Chooz started a search for the non-zero neutrino mixing angle #13 using the new gen- 
eration of v e detector since April 2011. A ratio of observed-to-predicted events of 0.944 ± 
0.016 (stat) ± 0.040 (syst) was observed in 101 days of detector running. We found sin 2 2#i3 
= 0.086 ± 0.041 (stat) ± 0.030 (syst) at Am^ m = 2.4 x 10" 3 eV 2 , based on rate and spectral 
shape information. The no oscillation hypothesis is excluded with 94.6% C.L., which can be 
interpreted as an indication of the non-zero #13. 



Table 1: Summary of uncertainties 



Source 


Uncertainty w.r.t signal 


Statistics 


1.6% 




Bugey4 measurement 


1/10/ 

1.4% 






.buel composition 


0.9 % 






Thermal power 


0.0 70 




Reactor 


Reference spectra 


U.O /o 


l.o /o 




Energy per fission 


U.Z /o 






IBD cross section 


U.Z 7o 






Baseline 


0.2 % 






Energy response 


L.I /0 






^delay containment 


0.6% 






Gd fraction 


0.6% 




Detector 


At e + n 


0.5% 


2.1% 




Spill in/out 


0.4% 






Trigger efficiency 


0.4% 






Target H 


0.3% 






Accidental 


< 0.1 % 




Background 


Fast neutron 


0.9% 


3.0% 




9 Li 


2.8% 
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Figure 3: Top: Observed and predicted positron energy 
pectation with stacked background spectrum. Bottom: 
difference between the best-fit sin 2 2#i3 



spectra for the best-fit sin 2 2613 and no-oscillation ex- 
Difference between data and no-oscillation (dots) and 
and no-oscillation expectation (line). 




Figure 4: Combined analysis with Double Chooz, T2K and MINOS assuming normal mass ordering. The A\ 2 
distribution for sin 2#i3 (Upper plot) and the allowed region in 8 vs sin 2613 plane (Lower plot), respectively. 



